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ABSTRACT 

We use the KOI- 13 transiting star-planet system as a test case for the recently 
developed BEER algorithm (Faigler & Mazeh 2011), aimed at identifying non- 
transiting low-mass companions by detecting the photometric variability induced 
by the companion along its orbit. Such photometric variability is generated by 
three mechanisms, including the beaming effect, tidal ellipsoidal distortion, and 
reflection/heating. We use data from three i^^ep^er quarters, from the first year of 
the mission, while ignoring measurements within the transit and occultation, and 
show that the planet's ephemeris is clearly detected. We fit for the amplitude 
of each of the three effects and use the beaming effect amplitude to estimate 
the planet's minimum mass, which results in MpSini — 9.2 ± 1.1 Mj (assum- 
ing the host star parameters derived by Szabo et al. 2011). Our results show 
that non-transiting star-planet systems similar to KOI-13.01 can be detected in 
Kepler data, including a measurement of the orbital ephemeris and the planet's 
minimum mass. Moreover, we derive a realistic estimate of the amplitudes un- 
certainties, and use it to show that data obtained during the entire lifetime of 
the Kepler mission, of 3.5 years, will allow detecting non-transiting close-in low- 
mass companions orbiting bright stars, down to the few Jupiter mass level. Data 
from the Kepler Extended Mission, if funded by NASA, will further improve the 
detection capabilities. 

Subject headings: methods: data analysis, planetary systems: detection, stars: 
individual: KOI-13 
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Introduction 



Traditionally, before the age of continuous high-precision space-based photometry, 
the investigation of low-mass transiting companions used spectroscopy primarily to study 
the system's orbit via the spectroscopic Doppler effect, and photometry primarily to 
study the system's geometric characteristics through measurements of eclipse (transit and 
occultation) light curve. However, it has been shown that spa ce-based phot ometry from 



missions such as Kepler fjBorucki et al. 



1998 



Auvergne et al. 



2010 



Koch et al. 



20101 ) and CoRoT fiRouan et al. 



20091 ) can be used to detect optical variability induc ed by a low-mass 



companion along its or 



2003 



pital motion, down to th e planetary mass level (e.g. 



Mazeh fc Faigler 



2010 



Welsh et al 



Jenkins &: Doyle 



20101). Hereafter we refer to the measurement of 



photometric variability correlated with the orbit as orbital photometry. 



Recently it was shown by 



Faigler fc Mazehl ( 1201 ll ) that analysis of orbital photometry 



can be used for detecting new non-eclipsing binary systems, where the secondary (star or 
planet) can be of low-mass, down to the planet mass range. Their BEER algorithm is 
aimed at identifying the system's ephemeris, and estimating the secondary minimum mass, 
M2sini, where M2 is the mass of the secondary an d i the angle betw e en th e orbital angular 
momentum axis and the line of sight. In addition, 
detection of two targets showing an orbital photometry signal, which in both cases suggests 
the presence of a non-eclipsing companion whose minimum mass is close to 70 Mj. We 
use here the same method, with some modifications an d additions, and app ly it to the 



Faigler fc Mazehl ( 1201 ll ) presented the 



recently discovered transiting planet system KOI-13.01 (IRowe et al. 



2011al lbl). containing 



a massive planet orbiting a bright A-type host star {Kp = 9.96 mag; KIC ID 9941662) 
ev ery 1.76 days. Since tha t system is transiting we use it as a test case for the approach 



of 



Faigler fc Mazehl (120111 ). We use only the out-of-eclipse light curve to show whether 



non-transiting star-planet systems similar to KOI-13.01 can be discovered in Kepler data. 
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We briefly describe the various mechanisms inducing photometric variabihty along the 
orbit in § [2J In § E] we present our data analysis where we show that the out-of-eclipse data 
can be used to clearly detect the ephemeris of KOI-13.01. In our discussion in §|l]we give 
an estimate of the companion minimum mass (which in this case is close to the true mass), 
and show that data from the entire 3.5 year Kepler mission lifetime is expected to allow 
detecting non-transiting close-in star-planet systems down to the Jupiter mass level. Our 
work is summarized in § [5l 



2. Orbital Photometry 



As described by iFaigler fc Mazehl (l201lh photometric variability along the orbit can 
be induced by three mechanisms: (I) beaming (aka Doppler boosting), (II) tidal ellipsoidal 
distortion, and (III) reflection of light from the star by the planet. The third mechanism 
includes the effect of heating of the planetary atmosphere by the star's radiation. 



The tidal ellipsoidal and reflection effects are well known from observations of eclipsing 



binaries, while the 



3 eamin g effect was observed o nly re cently. The 



Shakura fc Postnovl (119871) . 



Carter et al 



2011 



Maxted et al 



Loeb fc Gaudil f 2nn3h. and 



from both space (yan Kerkwiik et a. 



2000h . 



Faigler fc Mazeh 



20 







Zucke 



Mazeh fc Faiglei 



20 111 ) and the ground (IShporer et al 



r et al 



201G 



atter was described in 



(|2007h. and observe d 



Bloemen et al. 



2011 



2010: see also 



The equ ations describing t 



are given by 



le am plitudes, in relative flux, of each of the three effects 



Faigler fc Mazehl (120111 ). We list them here for completeness, for a Jupiter 
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mass companion: 

A.ea. = «bea.4^ = 2.7 ^— J J (^^^ J ppm , (1, 

M2 sin i f . . ^ „ . ■ f Rs\^ f Ms \ ^ f Port \ ^ f ^2 sin i 



Acmp = «emp^^ J sin . = 13 Ceiup sin ^ J J [j^ ) ) PPm 



(2) 



A,efl = a.efl0.1^:-iJ sinz = 57a,efisinz^gj (^1^ j (^^ J ppm, (3) 

where Ms and i?^ are the host star's (i.e. the primary star's) mass and radius, respectively, 
M2 and i?2 are the secondary mass and radius, respectively, Porb, i and a are the orbital 
period, inclination angle, and semimajor axis, respectively, K^y is the host star's radial 
velocity (RV) modulation semi amplitude, and c is the speed of light . The three coefficients: 



«beam, ^eiiip, and ctrefi are each of order unity — see 



Faigler fc Mazehl (120111 ) for more details. 



The expressions above are approximations, where the primary assumptions are a 
circular orbit, negligible luminosity for the low-mass companion relative to the primary, and 
M2 < Ms. 

Using the above three equations we show in Fig. [1] the expected amplitudes of the 
three effects for a star-planet system similar to KOI-13.01, and a range of orbital periods. 
The KOI-13.01 orbital period is marked by a vertical dashed line, where the expected 
amplitudes are an order of 10^^ to 10^*^ in relative flux. We note that Fig. [T]is plotted while 
assuming the light from the target is not blended on Kepler CCDs with light coming from 
nearby stars; such contamination will dilute the observed amplitudes. 



3. Data Analysis and Results 



We used here data from three Kepler quarters, Q2, Q3, and Q5. Data from Ql and 
Q4 were avoided as they were produced using insufficient aperture size to collect all the 
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flux of tliis saturated star, and data from QO (and Ql) sliow increased noise level relative 
to other quarters. We verified that using all data through Q5 gave identical results, albeit 
with larger error bars by a factor of ~ 2. 

Before applying our analysis we applied two preprocessing steps. First, we ignored the 
in-eclipse points from all transits and occultations, removing 18.3 % of the data. 

The raw Kepler data includes instrumental trends, temporal gaps and discontinuities 
(or "jumps") in the flux level. Correcting for these effects is important for detecting 
variability along the entire orbital phase, as opposed to the detection of eclipses whose 
typical duration is only a small fraction of the entire orbit. For that end we applied a second 
preprocessing step where we divided the Kepler light curve into several continuous sections 
without gaps or jumps. The typical duration of each such section is 1-2 months. We fitted 
each section with a 5th degree polynomial and divided by the fit to get the corrected light 
curve, in relative flux. The fitting was done iteratively while ignoring 5 cr outliers relative to 
the fit until no outliers were found. We checked that using a polynomial degree of 4 and 6 
gave consistent results, within 0.5 cr. This step removed an additional 5.6 % of the original 
three quarters data set, most of which from the start and end of the continuous sections. 

The corrected light curve is shown in Fig. [2] top panel. Data from QO, Ql and Q4, 
that were not used in our analysis, are shown in gray. As can be seen in the figure the data 
from QO and Ql, containing the first 45 days of data have a larger scatter than the rest. 
The data we use in the following analysis, from Q2, Q3, and Q5, are shown in black. They 
contain 9,689 Long Cadence measurements, each with an effective exposure time of 29.4 
minutes, spanning 368.4 days. Fig. [2] bottom planet shows a zoomed-in view of 10 days of 
data from Q5, where the orbital photometric variability can be visually identified. 

To look for the orbital photometric signal, including the identification of the orbital 
period and phase, and the measurement of the amplitudes of each of the photometric 



- 7- 



effects, we followed the 2-step BEER algorithm presented by iFaigler fc Mazehl (120111 ) 



Step 1: The goal of the first step is to detect the orbital period. To do so we constructed 



a periodogram where for each trial period we fitted, in a linear least 



Shporer fc Mazeh 



double harmonic model (e.g., 
parameters, or coefficients: 



/'2vr A . /27r , 

f{t) = ao + aic cos ( —t j + sm I —t 



2006 



Shporer et al. 



squa res manner, a 



20101 ) including 5 



a2c cos 



2n 



t' ] + a2s sin 



2n 



t' , (4) 



P/2 J \P/2 

where / is relative flux, P is a trial orbital period, and t' = t — To where t is time and Tq 
is the time zero point determining the phase. Tq is arbitrarily chosen in the first step and 
set to time of inferior conjunctioiij in the second step (see below). For each trial period the 
statistics we used is based on the statistics and is taken to be: 

Ax' 



Amean A 



(5) 



X X 

where Xmean is the for fitting a single parameter model, the measurements' mean. This 
method does not depend on an accurate knowledge of the measurement errors scaling. Our 
periodogram is shown in Fig. [3l The strongest peak corresponds to an orbital period of 
Porb = 1.7637 ± 0.0013 day. Peaks marked by a vertical solid black line are associated with 
Porb, including harmonics and a sub-harmonic. 

Step 2: Since the phase we used for the detection of the orbital period in the first step 
was arbitrary, the fitted coefficients do not necessarily have a one to one correspondence 
with the three physical effects described above. Therefore, in the second step we shifted Tq 
to have phase zero positioned at inferior conjunction. This was done by using the fitted 
coefficients for the best fit period obtained in Step 1, and requiring that (I) a2s, the only 
coefficient among the four that does not correspond to any physical effect, will vanish. 



^In a transiting system this is mid transit time. 
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and (II) a2c wil 
ambiguous (see 



be negative anc 




Faider & Mazeh 


2011) 



positive. The first requirement alone will still be 



20111). The derived transit time is Tq = 245 55138.7439 ± 



Borucki et al. 



0.0013 (BJD). The ephemeris derived here is consistent with the ephemeris of[ 
(120111) for this objec t, and they are both listed in Table [H The higher precision of the 



Borucki et al 



(l201l[ l ephemeris is to be expected as it is derived from data covering the 



entire orbital phase, including transit and occultation. 

After applying the shift we refitted the model using the period found in the first 
step, and the resulting fitted coefficients are listed in Table [2J A quick examination of the 
coefficients shows their signs are each as expected, meaning negative for the refiection (aic) 
and ellipsoidal (a2c) effects and positive for the beaming effect (ais). 

We tried to estimate realistic errors on the fitted coefficients. The linear least squares 
method provides us "formal" errors, but they do not account for possible remaining 
systematic trends in the Kepler light curve. To account for systematic trends we repeated 
our analysis while ap plying cyclic permutations to the residuals of our best fit model (e.g.. 



Carter fc Winn 



20091). Each of those permutations was generated by repeatedly shifting the 
residuals of every measurement and adding it to the model value of the subsequent point 
in the time series. The residual of the last measurement is added to the first point. Cyclic 
permutations keep the order of the residuals unchanged, preserving possible systematic 
trends, unlike random permutations. After analyzing all possible 10"^ cyclic permutations we 
examined the distributions of the resulting coefficients and took half the difference between 
the 84.13 and 15.87 percentiles to be the 1 a uncertainty. The coefficient's distributions 
were highly symmetric about the median, and the latter were indistinguishable from the 
original fitted values. The errors derived this way were close to twice as large as the errors 
obtained directly from the linear least squares method. 



Our phase folded light curve, based on the ephemeris derived here is presented in the 
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top panel of Fig. HI That figure does not include the in-eclipse data, emphasizing that it 
was not part of the analysis. The sinusoidal modulation constructed from the coefficients 
fitted here is overplotted with a solid line and the filled circles represent the binned light 
curve. Fig. H] bottom panel shows the individual sinusoidal effects plotted using the fitted 
coefficients, including the beaming (red), ellipsoidal (blue) and refiection (green) effects. 
The solid line (black) is the fitted model, the sum of the three modulations. The bottom 
panel shows that the overall double peak shape of the folded light curve results from the 
ellipsoidal effect, the large difference between the two minima results from the reflection 
effect and the small difference between the two maxima is the signature of the beaming 
effect. This panel shows clearly that the beaming effect amplitude is much smaller than 
the other two, which is expected according to Fig. [H although it is detected beyond a 10 cr 
signiflcance (see Table [2]). 



4. Discussion 

Dilution of the light from the target with light from a nearby star, or stars, will 
act to decrease the observed amplitudes of the three effects. Therefore, before using the 
fltted coefficients to extract physical information about the star-planet system one needs 
to account for possible dilution of light. For KOI-13 the dilution is significant since the 
target is the brighter member of a visual binary, with a separation of ~ 1 arcsec. Since 
Kepler's pixels span 3.98 arcsec the light from the two stars is completely blended together. 
To translate our measured photometric amplitudes into the undilut ed amplitude s , and in 



Szabo et al. 



(12011 



see 



turn measure the systems physical parameters, we use the results of 
their Table 1). We take their V band magnitude difference, of 0.29 mag, to calculate the 
dilution, while adopting an error of 0.1 mag, since no error is given. We use the V band 
magnitude difference as an approximated magnitude difference in the Kepler band, since 
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the two stars have similar T^fj, and they are located at the same distance and position in 
the sky. Therefore, the dilution factor, V, by which the coefficients need to be multiplied in 
order to derive the undiluted amplitudes is 1.77 ± 0.07. 

The amplitude of the beaming effect can be used to estimate the minimum mass of the 
secondary (M2sini). Rewriting Eq.[T] gives: 



Mo sin i 



0.37 / iVf, \ ^ f Porb\ ^ ( ^beam 



M.J . 



(6) 



We calculate the right hand side of Eq. E] using Port derived here, an undiluted beaming 
amplitude o f Ahpnm = l^^is = 9.32 ± 0.86 ppm, and a host star mass of Mg = 2.05 Mq from 



Szabo et al. 



( 1201 ll ). We adopt a 10 % uncertainty on the stellar mass since no error is given. 
The beaming coefficient is estimated to be abeam = 0.73 ± 0.03 assum ing the host star 



201ll ). and by integrating 



radiates as a blackbody with Tf.ff = 8,510 ± 390 K (iSzabo et al. 
over Kepler's transmission spectrum. The resulting minimum mass for the companion is 
M2sini = 9.2 ±1.1 Mj. We note that a different estimate of Mg and T^ff will result in 
a different M2sinz. 

Since this system is transiting the above minimum mass is close to the true mass, and 
perhaps consistent with the latter within the error bars. This puts the companion's mass 
at the massive planet mass range. 



As noted by 



Faider fc Mazehl (12011 



Kilic et al 



201l[ ). given sufficient 



see also 

information about the host star Eq. |5] shows that v4e//ip can be used to estimate another 
expression for the minimum mass, M2sin^i, and combined with M2sini estimated above 
the orbital inclination angle can be derived, which in turn allows to derive M2. Meaning 
that analysis of orbital photometry signals can, in principle, result in determining the 
companion's true mass and not only the minimum mass. However, while Eq. [H describing 



the amplitude of the beaming eff ect is fa.irly ac curate, Eq. |2l desc ribing the e 



amplitude is approximated (e.g.. 



Morris 



19851 ). It was shown by 



Pfahl et al 



lipsoidal effect 



(120081 ) that 
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the amplitude of this effect for an early type star with no convective layer, like KOI- 13, is 
difficult to describe analytically and it is likely far from the familiar expression used in Eq. [2] 
(although the shape of the effect can still be described as a cosine at the first harmonic of 
the orbital period). Therefore, we do not proceed here with using Aeiup to derive M2 and 
i. The lack of an adequate equation for AeUip does not allow us to use the ratio of the 
elhpsoidal and beaming amplitudes to check the likehhood of the model, as described by 



Faider fc Mazehl (1201 ih . 



The error bars derived for the fitted coefficients (see Table [2]) show that the KOI-13 
light curve is sensitive to modulations down to 2 ppm which can be detected at the 3-4 a 
level. Assuming white noise, the data from the entire 3.5 years Kepler mission lifetime can 
be sensitive to 1 ppm modulations for a target similar to KOI-13 in brightness and stellar 
noise. It is interesting to use these sensitivities to calculate the companion's mass that can 
be identified for a range of orbital periods. To derive this predicted sensitivity of a search 
for orbital photometry signals in Kepler data we adopt a signal to noise threshold which is a 
factor of two larger than the above, bringing it to 7 a. Although our error bars were already 
inflated by a factor of two relative to the formal error bars (see §[3]), this additional factor 
of two is required in order to make th e amount of fa lse positives small enough considering 



the large number of light curves (e.g. 



Jenkins 



2OO2I ). Such sensitivity diagrams are shown 



in Fig. [5l The left panel is for a KOI-13-like host, and the right panel is for a Sun-like host 
at Kp=12 mag. Only the beaming and ellipsoidal effects are plotted on those diagrams, 
since they are sensitive to the companion's mass. The reflection effect is sensitive to the 
companion's radius (see Eq. [3]), which is only weakly dependent on mass at the mass range 
of gas giant planets and brown dwarfs. Both diagrams assume no dilution and an arbitrarily 
chosen orbital inclination angle of i = 60 deg. The lines in both panels correspond to the 
expected sensitivity once data from the entire mission lifetime is obtained, taken to be 2 
ppm in the left panel and 5 ppm in the right panel, detected at the 7 a level. 
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A close inspection of Fig. [5] shows that orbital photometry is sensitive to short-period 
low-mass companions down to the few Jupiter mass level. Therefore, it can be used to 
detect low-mass stellar companions, brown dwarfs and massive planets. This mass range 
completely overlaps the so called bro wn dwarf desert — a paucity in the mass function 



of companions to Sun-like stars (e.g.. 



Grether fc Lineweaver 



2006 



Sahlmann et al. 



2011 



Evans et al. 



20111 ). Occupants of the brown dwarf desert are rare, but orbital photometry 



can detect most if not all short-period brown dwarfs orbiting bright stars in the Kepler field, 
thereby characterizing the high mass end of the planetary mass distribution and the low 
mass end of stellar companions mass distribution. 

The assumption of white noise above and the sensitivities presented in Fig. |5] could 
be somewhat optimistic. Later type stars, colder than KOI-13, constitute the majority 
of Kepler targets and are expected to show activity at a level that may exceed the 
Kepler precision. This will degrade the sensitivity of the BEER algorithm to the orbital 
photometric signal. Stellar oscillations also represent a potential problem for the BEER 
approach, although unless the oscillating frequency is closely related to the orbital frequency 
it is possible they can both be resolved. Both stellar activity and oscillations could 
in principle mimic an orbital photometric signal, resulting in a false BEER detection. 
Therefore, BEER detections should be examined ca refully, specifically that the coefficients' 



Faigler fc Mazeh 



201l|). In addition. 



sign and ratio are as expected (for more details see 
many of the Kepler targets are blended with other nearby stars on Kepler CCDs, as 
is KOI-13 itself. Blending will decrease the overall amplitude of the orbital signal but 
otherwise will not affect its shape, resulting in underestimating the companion's mass. All 
the above emphasize the need for RV follow-up to confirm BEER detections and measure 
the companion's mass. 



On a more positive note, the errors listed in Table |2] were obtained from an incomplete 
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light curve, without the transit and occultation phases which cover 18 % of KOI-13.01 
orbit. Therefore, the results obtained here and the sensitivities shown in Fig. E] should be 
somewhat better for a non-transiting system similar to KOl-13.01, where the entire light 
curve is used, unless the system's orientation is close to face-on. 

Analysis of orbital photometry can, in principle, be used to look for additional signals, 
induced by additional objects in the system, including a possible second planet which is not 
necessarily transiting. This possibility is intriguin g in light of the large number of candidate 



multi transiting planet systems found by Kepler fiLissauer et al.l 120111 ) . and the fact that 



about one sixth o f the candidate tran sit ing planetary s ystern s contain more than a single 



transiting planet ( iBorucki et al. 



201ll ). iLatham et al.l (120111 ) have shown that massive 



planets, like KOI-13.01, are less likely to reside in multi transiting planet systems, although 
it is still possible that additional planets occupy non-transiting orbits. 

Our periodogram presented in Fig. [3] shows two small peaks, marked by dashed vertical 
lines, that are not associated with the 1.76 day orbital period. The strongest among the 
two corresponds to a period of 1.0600 ± 0.0010 days, and the smaller peak is at the first 
sub-harmonic. We verified this secondary photometric signal, whose semi amplitude is 12 
ppm, exists in each of the three Kepler quarters included in our analysis. 

We investigated this further by subtracting the 1.76 day double harmonic model and 
applying the same period analysis as described above to the residuals. The derived fitted 
coefficients were af^ = 2.3 ± 1.0 ppm, a^^J = 11.8 ± 1.0 ppm, and a'^^ = —3.1 ± 0.5 ppm, 

(2) 

where the upper index is used for marking a secondary photometric signal. The sign of a\J 
and the relative size of the coefficients makes it unlikely this signal originates from a second 
planet orbiting at 1.06 days. 

The ratio between the 1.06 day period and the orbital period is 0.60100 ± 0.00072, 
which is close to an integer numbers ratio of 3:5. Assuming this is not a coincidence, it 
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raises the possibility that the 1.06 day signal is a stellar puls ation triggered by the p 



anet ary 



2OIOI), a 



companion. A similar scenario was suggested for WASP-33 (iCoUier Cameron et al.l 
planetary system that resembles KOI-13 as it includes a short period (1.22 days) transiting 
planet orbiti ng a pulsating A-type star whose spin axis is rn isaligned with the planetary 



orbital axis (jCoUier Cameron et al. 



2010 



Barnes et al 



20111 ). An extr eme example of 



2Qll|), 



pulsations induced by a binary companion is the KOI-54 binary system (1 Welsh et al. 
composed of two similar A-type stars on a highly eccentric orbit. Still, this scenario does 
not explain why we see a periodicity specifically at a relation of 3:5 with the orbit and not 
another integer relation. 

We conclude that the true nature of this secondary photometric signal is currently not 
clear, although unlikely to be due to a second planet. In addition, we note that since the 
target is blended with a nearby star on Kepler CCDs, this secondary photometric signal 
could originate from the nearby star. A more detailed study of this signal, once more 
Kepler data becomes available, will be important for better understanding the possible false 
alarm scenarios of the application of the BEER approach to Kepler data, which have not 
been thoroughly studied yet. 

Finally, we note that a close inspection of the phase folded light curve in Fig. H] top 
panel shows there are small deviations of the binned light curve from our double harmonic 
model. That model is simplistic, meant to be used only as a detection tool, and not as a 
detailed physical model, especially for data of such high quality. Therefore, the presence of 
these deviations should not be alarming, but rather be viewed as an opportunity for a more 
detailed study of bright systems using Kepler data. 
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5. Summary 



We have pre sented here a test case for analyzing orbital photometry using the 



BEER approach (jPaigler fc Mazeh 



20111 ). We have shown that photometric measurements 
taken out of eclipse (transit and occultation) can be used to detect the presence of 
KOI-13.01, measure its orbital ephemeris and estimate its minimum mass. Therefore, 
similar non-transiting systems can be detected using Kepler data. 

Moreover, we used the error bars derived for the fitted coefficients and a 7 cr detection 
threshold to estimate the sensitivity of Kepler data for low-mass non-transiting companions. 
Assuming white noise we showed that data accumulated during the entire 3.5 years 
Kepler mission lifetime will allow detection of low-mass short-period companions down to 
the few Jupiter mass level, for bright stars. In case the Kepler mission will be extended, 
this sensitivity will improve, allowing detection of lower mass companions, orbiting fainter 
stars at longer periods. 

Spectra of hot stars, of early spectral type like KOI- 13, are not rich with deep and 
narrow absorption lines which are used for measuring high precision RVs. As a result 
those stars are not regularly monitored by RV surveys, leading in turn to highly limited 
knowledge about their low-mass companion population. Analysis of the orbital photometry 
signal presents an opportunity to extend our knowledge about low-mass companions to 
early type host stars, especially considering they show a decreased amoun t of stellar acti vity. 
For ex ample, orbital photometry can be used to test the suggestion of 



Bouchy et al. 



( 120111 ). that massive planets and brown dwarfs are more frequent around F-type stars 
than G-type stars, and see if this trend continues into A-type stars. RV confirmation of 
such companions, detected through orbital photometry, will probably be difficult or even 



impossible (For an RV confirmation o: 



photometry see 



Ehrenreich et al. 



a white dwarf companion detected through orbital 



20111). However, a detection of a large sample along with 
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an improved understanding of the false positive scenarios will allow a statistical analysis 
and a comparison to the characteristics of the companion population orbiting Sun-like stars. 

The analysis done here uses only out of eclipse data. Several authors have already 
shown that Kepler's continuous and high precision photometr y taken in-eclipse can be 



used for additional characterization of the companion's o rbit ( iBarnes 



2009 



Barnes et al. 



2011 



Shporer et al 



2011 



Groot 



2011 



Szabo et al. 



20111 ). These authors show that 



small asymmetries in the eclipse light curve can be used for measuring the sky-projected 
spin-orbit angle — angle between the star's rotation axis and the companion's orbital 
angular momentum axis. Thereby showing the high scientific potential of using photometry, 
out and in eclipse, to study binary systems, those with low-mass companions in particular. 

Kepler was competitively selected as the tenth Discovery mission. Funding for this 
mission is provided by NASAs Science Mission Directorate. We warmly thank Michael Endl 
for useful comments. A.S. acknowledges support from NASA Grant Number NNX10AG02A. 
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Table 1. KOI-13.01 ephemeris 



Parameter This work Borucki et al. (2011_) 



Orbital period, Porb (days) 1.7637 ± 0.0013 1.7635892 ± 0.0000014 

Inferior conjunction time, Tq (BJD) 2455138.7439 ± 0.0013 2454953.56498 ± 0.00012 
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Table 2. Fitted coefficients 



Coefficient Effect Value 

[ppm] 



aic Reffection -39.78 ± 0.52 

ais Beaming 5.28 ± 0.44 

a2c Effipsoidal -30.25 ± 0.62 

a2s — 0.0 ± 0.48 




Period [day] 



Fig. 1. — Expected amplitudes of the beaming (red), ellipsoidal (blue), and reflection (green) 
effects for a range of orbital periods, for a system similar to the KOI-13.01 star planet system. 
The plot is in log-log scale and is based on Eqs. [THSl The dashed black line marks KOI-13.01 
orbital period. 
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Fig. 2. — Top: KOI-13 Kepler Long Cadence light curve, after applying the steps described 
in §[3l including removing the transit and occultation data. The data from QO, Ql and 
Q4 are shown in gray and were not used here. The data from Q2, Q3 and Q5 are shown 
in black. Bottom: Zoomed-in view of a 10 day light curve during Q5 (without in-eclipse 
measurements) . 
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Fig. 3. — The double harmonic periodogram (blue), where the strongest peak corresponds 
to the orbital period, of 1.76 days. The solid vertical lines (black) at the top of the figure 
mark frequencies which are a sub harmonic or harmonics of the strongest peak. The dashed 
vertical lines (black) mark unidentified frequencies where one is an harmonic of the other. 
See text for more details. 
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Fig. 4. — Top: phase folded light curve, using the ephemeris derived here. The gray dots are 
Kepler Long Cadence data and the black circles are the binned light curve (error bars are 
smaller than the markers). The solid line represents the sinusoidal modulation constructed 
using the fitted coefficients. Bottom: we show the sinusoidal signals, in dashed lines, of each 
of the three effects plotted using the amplitudes found here (B=Beaming, E=Ellipsoidal, 
R=Reflection). The solid line (black) is the fitted double harmonic model, which is the sum 
of the three effects. 
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Fig. 5. — Minimum detectable secondary (star, brown dwarf or planet) mass as a function of 
orbital period for a 7 o" detection, and assuming an inclination angle of 60 deg. The reflection 
effect can not be put on these diagrams, since its amplitude is sensitive to the companion's 
radius, which in turn is not sensitive to the companion's mass in the mass range of gas giants 
and brown dwarfs. Left: Assuming a KOI-13-like host (primary) star, and an amplitude of 
2 ppm detected with Kepler data from the entire 3.5 years mission lifetime. Right: Similar 
to the left panel, but for a Sun-like host star, at a brightness of Kp = 12 mag and assuming 
5 ppm amphtude. 
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